RBRC 941 



KNO scaling of fluctuations in pp and pA, and eccentricities in heavy-ion collisions 

Adrian Dumitru 1, 2 and Yasushi Nara 3 

1 RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY- 11 973, USA 
2 Department of Natural Sciences, Baruch College, CUNY, 
17 Lexington Avenue, New York, NY 10010, USA 
3 Akita International University, Yuwa, Akita-city 010-1292, Japan 

Multiplicity fluctuations at midrapidity in pp collisions at high energies are described by a negative 
binomial distribution and exhibit approximate Koba-Nielsen-Olesen (KNO) scaling. We find that 
these KNO fluctuations are important also for reproducing the multiplicity distribution in d + Au 
collisions observed at RHIC, adding to the Glauber fluctuations of the number of binary collisions 
or participants. We predict that the multiplicity distribution in p + Pb collisions at the LHC also 
deviates little from the KNO scaling function. Finally, we analyze various moments of the eccen- 
tricity of the collision zone in A + A collisions at RHIC and LHC and find that particle production 
fluctuations increase fluctuation dominated moments such as the triangularity £3 substantially. 

$h I Charged particle multiplicity distributions in the central region of inelastic (non-single diffractive) p+p collisions at 
high energies were shown by the UA1 and UA5 collaborations to follow a negative binomial distribution (NBD) p], Q 
exhibiting approximate "KNO scaling" Q over at least a limited range of multiplicities excluding the tails. Within 
the framework of high-energy QCD they may be thought to arise from fluctuations of the density of large- a; valence 
charges [3| and of stochastic emissions in the rapidity evolution ladders 0, @ leading from the rapidity of the sources 
to the central region. 

Collisions of hadrons or heavy ions at high energies release a large number of gluons from their wave functions. In 
fact, the wave function of a hadron boosted to (nearly) the light cone is so densely packed with gluons that they may 
"overlap", leading to non- linear interactions Q- Therefore, at high energies the colliding hadrons can be treated as a 
high occupancy gluon field. This dense system is nowadays referred to as Color Glass Condensate (CGC) Q. 
Here, we use the "fc_L-factorization" approach Q to compute particle production in high-energy collisions: 
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where N c — 3 is the number of colors in QCD and K ~ 1.5 — 2 is a multiplicative factor to account for corrections 
to this LO formula. Further, we assume that the hadron multiplicity is proportional to the multiplicity of gluons. 
These factors may depend somewhat on the initial condition for small-s evolution of $(fcj_) but were found to be 
approximately independent of collision energy or centrality (for heavy-ion collisions) Q. Also, in eq. (fTJ) £1,2 = 
(p±/y/s) exp±y and the strong coupling is evaluated at the scale Q — max(|px + kx|, |p_L — kx|)/2. 

We require the evolution of the so-called unintegrated gluon distribution $>(k±,x) (per unit transverse area) with 
the light-cone momentum fraction x, starting from an initial condition at xq ~ 10~ 2 . This is obtained by solv- 
ing the non-linear Balitsky-Kovchegov (BK) equation [llj with the running-coupling kernel according to Balitsky's 
k> . prescription [±2j]. Specifically, we use the unintegrated gluon distribution "set MV" from ref. 0. 

For the case of heavy-ion projectiles and/or targets, we allow for fluctuations of the locations of the sources (i.e., 
of the valence charges at x$) for the small-x fields in the transverse plane before the collision [ij [l(| • This leads to 
fluctuations of the "geometry" of the collision zone from configuration to configuration, and to fluctuations of the 
number of participants iV par t and the number of collisions iVcoii which are determined within the well-known Glauber 
approach. Note that eq. (JTJ) refers to a single such configuration. We computed these "geometry" fluctuations assuming 
that the hard valence charges are smeared over a finite and energy independent area <7o ~ cr./v./v(200 GeV) = 4.2 fm 2 . 
This reduces higher-order eccentricities as compared to point-like nucleons [l3| which are used in some Monte-Carlo 
Glauber simulations. Our numerical simulations do not account for correlations (in the transverse plane) among the 
valence charges which could further suppress geometry fluctuations (l4l |. 

The unintegrated gluon densities Q(k±,x) from eq. |U have already been averaged over the local fluctuations of 
the valence charges in color space, and over the evolution ladders. It is in this sense that we interpret eq. (fljl as a 
mean (local) multiplicity. In each cell A 2 rx of the transverse plane the actual multiplicity is a NBD random variable, 

P(n) = F ( fc + ") ( 2) 

1 ; r(fc)r(n + 1) (n + k) n + k ' 1 ; 

Here, n = (dN/dr) d 2 r±) A 2 r±Ai] is the mean multiplicity from eq. (flj in a given cell and k is the fluctuation 
parameter: smaller k correspond to larger fluctuations about the mean and KNO scaling is obtained when k <C n 
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(see below). We finally average over geometric configurations of sources in the rj_ plane described above, and over 
the impact parameter of the collision. 

There have been numerous theoretical discussions of multiplicity fluctuations in high-energy collisions. Ref. [4j, in 
particular, argued that NBD multiplicity fluctuations arise in a semi-classical calculation of gluon production from 
dense valence charge sources. They obtain that the fluctuation parameter k is proportional to the density per unit 
transverse area of valence charge squared, i.e. to the saturation momentum Q 2 S at xq- 
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FIG. 1: (Color online) Left: multiplicity distribution of charged particles at \rj\ < 0.5 in pp collisions at yfs — 900 GeV. 
Stars show the result of our calculation (see text) while solid and dashed histograms correspond to data taken by the ALICE 
collaboration with the "NSD" and "INEL" triggers, respectively [l^]. The dashed vertical lines indicate the average and two 
times the average multiplicity, respectively. Center: Same at i/s = 2360 GeV. Right: Same at i/s = 7000 GeV compared to 
CMS NSD data EJl. 



We first analyze the multiplicity distributions in proton-proton collisions at LHC energies (fig. [TJ. We concentrate 
on the bulk of the distributions, AT c h < 3(iV c h) where (N c h) denotes the average charged particle multiplicity at a 



given energy. Over this range the data can be described reasonably well by a NBD with constant 
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Here At] = 1 and A 2 r^ is the area of a cell in transverse coordinate space over which we integrate eq. |TJ). Also, we 
choose Aqcd = 0.24 GeV. Numerically, kjn ~ 0.16 for p + p collisions at 2.36 TeV. We have checked that a weak 
energy dependence of k is allowed as long as it does not change the distribution P(A r c h) appreciably over the range 
that we are interested in. The tails of P(A c h) could be more sensitive to the detailed dependence of k on energy 1 but 
we do not explore the region A ch > 3(A ch ) here; see, for example, ref. [l7| . 

The most important consequence from @ is that since fc=const and smaller than the average multiplicity n, it 
follows that our multiplicity distributions satisfy Koba-Nielsen-Olesen (KNO) scaling Q. That is, the probability 
distribution P(N c h) is independent of energy if expressed in terms of z = AT c h/(iV c h); for n 3> k and in the region 
z> k/n the NBD © can be written in the form of a Gamma distribution 
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We show the KNO scaling function in pp collisions explicitly in fig. [3] below. 

In fig.[5]we compare the calculated charged particle multiplicity distribution in d+Au collisions at </s — 200 GeV to 
uncorrected data from STAR [ijj]. As described above, here we include also fluctuations of the number of participants 
Afpart and of the number of binary collisions N co n which arise for different configurations of nucleons in the target 
nucleus. Within our formalism, N co \\ fluctuations alone are insufficient to reproduce the experimental multiplicity 
distribution. In this case we obtain a peak in P(A r c j 1 ) before the cutoff of the distribution which can be traced back 
to the fact that A^ c h does not increase linearly with the density of sources when the latter is high. This "saturation" 
of particle production is also responsible for the higher elliptic eccentricity of the collision zone than obtained from 
simple linear estimates (l9j . 

Additional intrinsic fluctuations with 
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1 The same applies to rapidity intervals bigger than < 0.5; UA5 found that k then actually decreases with energy 0. 
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FIG. 2: (Color online) Left: multiplicity distribution of charged particles at < 0.5 in min. bias d + Au collisions at 
= 200 GeV. Various models (see text) are compared to uncorrected data from the STAR collaboration (circles) [T^ ]. 
Right: multiplicity distribution predicted for min. bias p + Pb collisions at y/s = 4400 GeV. 



lead to a good fit to the data; such scaling of k with the number of sources is expected due to the way that negative 
binomial distributions add 2 . On the other hand, fcd+Au ~ max (Ta(t±), Tg(rx)) produces a multiplicity distribution 
inbetween the above cases, exhibiting too little fluctuations. Once again, it is reasonable that the magnitude of 
fluctuations is determined mostly by the dilute source (as also assumed in ref [13 )• Our prediction for p + Pb 
collisions at LHC is shown in fig. [2]on the right; this corresponds to (N^) ~ 16 and k from eq. (A prediction for 
the multiplicity distribution in p + Pb collisions at the LHC from the "KLN model" was shown previously in ref. 20]). 
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FIG. 3: (Color online) KNO scaling plot of the multiplicity distributions of charged particles at 77] < 0.5. 

Due to the presence of iVcoll fluctuations our multiplicity distribution for p + Pb does not exhibit exact KNO scaling, 
as seen in fig. [3] Nevertheless, for \rj\ < 0.5 and N c h < 3(A c h) we predict relatively small deviations from the KNO 
scaling function determined from p + p collisions. This is an important check for the presence of strong intrinsic 
particle production fluctuations (at fixed A part and A^ co ii) for a heavy-ion target. 

We now proceed to discuss the relevance of particle production fluctuations for various harmonic moments of the 
"eccentricity" of gluons produced in the initial state of heavy-ion collisions. We define moments of the initial density 



2 If x and y are two random variables with a negative binomial distribution with mean /i and fluctuation parameter k then z = x + y also 
follows a negative binomial distribution with mean 2fi and k z = 2k. Hence k is an extensive quantity proportional to volume, just as n. 
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distribution (preceding the hydrodynamic expansion in A + A collisions) in terms of the eccentricities [2lj, [22( 

(6) 



\J (r 2 cos nip) 2 + (r 2 sin nq 



(r 2 ) 

Other definitions are sometimes also used in the literature, see for example [23l - [25| . (•) denotes an average over the 
distribution of produced gluons in the transverse plane, dN/drid 2 r±- and rj_ = r(cos cp, sin <p). 

The eccentricities e„ are of interest because through hydrodynamic response they generate the flow harmonics and 
angular correlations in the final state of heavy-ion collisions [2l| - l33| . (Fluctuations in small- a; evolution may also 
lead to detectable azimuthal momentum anisotropies in high-multiplicity pp collisions at the LHC @ which are not 
due to "flow"). Flow harmonics in heavy- ion collisions have been published by the PHENIX [34[ and ALICE (35| 
collaborations. 




FIG. 4: (Color online) Centrality dependence of various moments e n of the eccentricity in Au + Au collisions at \fs = 
200 GeV. Dotted lines correspond to local fluctuations of iVcoll ("geometry fluctuations") only; dashed and full lines add 
particle production fluctuations according to a negative binomial distribution with k — k vv = const or k ~ min (Ta,Tb), 
respectively. 



In fig. [4] we compare the centrality dependence of £2 - £5 for three different models. In all cases, the lowest 
curve corresponds to the model with "geometry fluctuations" only, as usually considered in the literature. If particle 
production fluctuations according to a negative binomial distribution are added then in general e„ increases. The 
effect on £2 at N palt < 250 is small since the ellipticity of the overlap zone is of course dominated by the "almond 
shaped" geometry of heavy-ion collisions at finite impact parameter. On the other hand, we observe a large increase of 
e n for all n > 3 over the entire range of impact parameters, as well as of €2 for A^ part — > 2A, since these observables are 
fluctuation dominated. Most importantly, the ratio e^/e2 in mid-central collisions increases significantly. The largest 
increase is obtained if the fluctuation parameter k — k pp does not increase with the density of sources. Even for the 
more realistic case where k ~ min(TA,Te), higher-order eccentricities can increase by as much as 50%. We mention 
also that simulations using the DIPSY Monte-Carlo which performs the small-x dipole evolution stochastically have 
predicted a large €3 (36|, although the relation to KNO scaling in pp and pA collisions at the LHC had not been 
pointed out. 
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To summarize our main results: we found that in order to reproduce the measured multiplicity distribution in d+Au 
collisions at RHIC within the CGC approach it is important to take into account particle production fluctuations 
(according to a negative binomial distribution). We predict that these dominate over Glauber fluctuations also for 
p + Pb collisions at the LHC, resulting in a multiplicity distribution which is close to the KNO scaling function 
measured in p + p collisions. The effect of particle production fluctuations can be large also for some observables in 
heavy-ion collisions, such as for higher-order eccentricities. It will be interesting to see how this reflects in higher-order 
flow coefficients predicted by viscous hydrodynamics or in the centrality dependence of the jet quenching parameter 
Raa(p±) HI- 
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Appendix A: Eccentricities for heavy-ion collisions at LHC energies 
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FIG. 5: (Color online) Centrality dependence of various moments e n of the eccentricity in Pb + Pb collisions at ^/s = 2.76 TeV. 
Dotted lines correspond to local fluctuations of iV C oH ("geometry fluctuations") only; dashed and full lines add particle produc- 
tion fluctuations according to a negative binomial distribution with k = k pp — const or k ~ min(rA,Ts), respectively. 



